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Here in this work we propose a modified gravity with the action of f{R) = y R 2 
Einstein-Hilbert action to describe the late time acceleration of the universe. We obtain the equation 
of the modified gravity both in the metric and Palatini formalisms. An asymptotic solution of gravity 
equations corresponding to a constant Ricci scalar causes a late time acceleration of the universe. 
We do a conformal transformation in the action of the modified gravity and obtain the equivalent 
minimally coupled scalar-tensor gravity. The equivalent Brans-Dicke gravity of this model is also 
studied. To examine this model with the observation, the perihelion Precession of the Mercury is 
compared with our prediction and we put an upper constraint of Ro < Hq. This range for Ro is 
also in agreement with the cosmological acceleration at the present time. Finally we show that this 
action has instability for the small perturbations of the metric in vacuum solution in which adding 
a quadratic term of Ricci scalar can stabilize it. 



I. INTRODUCTION 

Recent cosmological data from the Supernova Type 
la (SNIa) and Cosmic Microwave Background (CMB)[l|, 
0, 0] indicate that universe is currently in the acceler- 
ating phase of expansion. The effect of acceleration of 
the universe on the Supernova Type la (SNIa) is making 
their apparent luminosity dimmer than our expectation 
[HH 0- The cosmological constant can explain this ob- 
servational effect. The best fit of the data with the stan- 
dard (Friedmann Robertson Walker) FRW model results 
an energy density of 10~ 47 Gev 4 for the cosmological con- 
stant. This value is a 55 order of magnitude smaller than 
the energy density of vacuum in quantum field theory 
(for review see]7, 8]). The cosmological constant is a ge- 
ometrical term, however it can be regarded as a perfect 
fluid with the equation of state of w = — 1. This equation 
of state provides a constant density during the expansion 
of the universe. To have the present universe the cosmo- 
logical constant should be fine-tuned and a small change 
in this term induces a dramatic change in the destiny of 
universe. 

To solve the cosmological constant problem, the 
Quintessence model as a time varying cosmological con- 
stant due to slow rolling of a scalar field is proposed. 
This model can cure the fine-tuning problem at the early 
universe, however there is no phys ical candidate for the 
origin of this scalar field @, M, El El 03 0, El El, El 
EllH HI HH ■ The other alternative approach dealing 
with the acceleration problem of the universe is chang- 
ing the gravity law through the modification of action of 
gravity by means of using f(R) instead of the Einstein- 
Hilbert action. Some of these models as 1/R and loga- 
rithmic models provide an acceleration for the universe 
at the present time [H, [H HJ [H [H [H H Hi] . 

The back-reaction of the structures at large scales on 
the dynamics of the background of the universe [3fjj is 
also proposed as a candidate for interpretation of the 
acceleration of the universe. Structured FRW universe 
is another approach which provides a larger luminosity 
distance in the expanding universe, making the lumi- 



nosity of SNIa dimmer [3l|. Here in this work we pro- 
pose a model for the modified gravity with the action of 
f{R) = (R 2 — Rq) 1 / 2 to provide a positive acceleration for 
the universe at the later times of the cosmic history. In 
section [II] we derive the modified Einstein equations both 
in the metric and Palatini formalisms and show that the 
solution of maximally symmetric static universe could be 
a de-Sitter space. In section Hill we show the equivalence 
of the modified gravity with the scalar tensor models such 
as the the Brans-Dicke and non minimally coupled con- 
formal transformed scalar tensor theories. In Section ITVl 
we show the instability of this action for small perturba- 
tions of curvature around the vacuum solution and show 
that adding an R 2 term can stabilize this model. The 
conclusions are presented in Section fVl 



II. FIELD EQUATION OF f(R) = yip 
MODIFIED GRAVITY 
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For an arbitrary action of the gravity there are two 
main approaches to extract the field equations. The first 
one is the so-called "metric formalism" in which the vari- 
ation of action is performed with respect to the metric. 
In the second approach, "Palatini formalism", the con- 
nection and metric are considered independent of each 
other and we have to do variation for those two param- 
eters independently. Here we obtain the field equations 
in the two mentioned approaches: 



A. metric formalism 

A generalized form of the action due to geometrical 
part of gravity can be written as follows: 



(1) 



where in the case of f(R) = R we have the ordinary 
Einstein-Hilbert gravity. Variation of the action with re- 
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spcct to the metric results in the field equations as: By substituting Ricci scalar and tensor in equation ([3]) 

_^ we will have a forth order non-linear differential equation 

}'{R)R l _ LV —-f(R)g l _ lv —(V ll V v ~g ll yV a V a )f'(R) = SnGT^. for A(r) with no analytical solution. A simple case is the 
2 ^ vacuum solution which results in: 

Here we c hoose the geometric part of the Lagrangian as 
f{R) = \/R 2 — Rq which implies the field equations as: 



RRuv 



VR 2 - Rl 



Rq 



(R 2 - Rl) 3 /< 



-^R 2 -R 2 g^ 



(3) 



3i? 



R 2 -R 2 



(V M i?V„i? - g^V a RV a R)} = 8irGT^. 



For the special case of Rq = we will recover the familiar 
Einstein field equations. One of the features of this action 
is that we will have an intrinsic minimum curvature for 
the space. Taking the trace of ((3]) results: 



R 2 



VR 2 - R'i 

9R 



2x1 R 2 - R 2 



Rf> 



{R 2 -R 2 y 



R 2 -R 2 M 



V^RV^R] = kT. 



-[3V a W a R 
(4) 



As a special solution we consider a maximally symmetric 
space-time without the energy momentum source where 
R is independent of space and time. This condition re- 
quires R = v2Rq and it can happen at the ultimate stage 
of the expansion of the universe. Now we solve the field 
equations for the modified gravity in the two special cases 
of static spherically symmetric and spatially homogenous 
spaces: 



1. Spherical static space 

We take a spherically symmetric Schwarzschild-like 
metric as: 



ds 1 = 



- x ^dt 2 + e x(r) dr 2 + r 2 {d9 2 + sin(6l) 2 cZ</> 2 ) (5) 



It is just a straightforward task to find the Christoffel 
symbols related to this metric. The non-vanishing com- 
ponents of Ricci tensor are as following: 



A"(r) A' 2 (r) A'(r) 



Rn — 



R 2 



- A «[l-rA'(r)] + l 



-^33 

sin 9 



P - P -X(r)( A» A» A» 



where the Ricci scalar obtain as: 

fl = e-*M(A"-A' 2 + — - 1) + % 



(6) 



(7) 



(8) 



(9) 



R = V2R (10) 
Substituting in equation ^ we obtain A(r) as: 

Kr) = -H^ + --- 2 -^Ror 2 ) (11) 
r r z 12 

In analogy with the Schwarzschild metric in the Newto- 
nian limit we set c\ = — 2G c ^ and C2 = and we have 
the metric as follows: 



ds z 



2GM V2 2 2 2 



err 
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2GM V2„ Xj 2 
err 12 



r 2 {d0 2 + sin d 2 d^2) 



This metric in the generalized gravity is similar to the 
Schwarzschild-de Sitter space in the Einstein-Hilbert ac- 
tion where i?o = 2\f2 A plays the role of the cosmological 
constant. A the solution of generalized Einstein equa- 
tion for a perturbation around the vacuum solution in 
the spherically symmetric metric is given in Appendix. 

There is a list of the observations from the solar sys- 
tem to the cosmological scales the modified gravity can 
be examined. Here we use the perihelion of the Mer- 
cury to check the behavior of the orbit in this model and 
put constraint on R as the parameter of the generalized 
gravity. Mercury is the inner most of the four terres- 
trial planets in Solar system, moving with high velocity 
in Sun's gravitational field. That is why mercury offers 
unique possibilities for testing general relativity and ex- 
ploring the limits of alternative theories of gravitation. 
The observed advance of the perihelion of Mercury that 
is unexplained by Newtonian planetary perturbations or 
solar oblately is [32] ]: 

Auj bs — 42.980 ± 0.002 arc-second/ century (13) 
= 2?r(7.98734± 0.00037) x IQ^radians/ revolution 

The calculation of the perihelion of Mercury in the 
Schwarzschild-de Sitter metric has been studied so far 
in [33| . For the time like geodesies they used Jacobi's 
inversion problem and found the constraint of A < 
10~ 55 cm~ 2 . From the spherical solution we have the sim- 
ilar constraint of Rq < 10 -55 cm~ 2 in our model. Com- 
paring with the horizon size of universe at the present 
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time implies i? < Hq 



2. cosmological solution 

We start the cosmological solution with the dynamics 
of universe in the radiation dominant epoch. For the 
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FIG. 1: Dynamics of y = R/Ro in terms of r = tR for the 
radiation dominate epoch. This dynamics results from the 
metric formalism for the field equation. The action for the 
modified gravity is taken as f(R) = y/ R 2 — Rq. 



does not alter the dynamics of the early universe from 
that of Einstein-Hilbcrt action. 



B. Palatini formalism 

In the Palatini formalism, the connections are con- 
sidered as independent variables, different than the 
Christoffel symbols of the metric [35[. Varying the ac- 
tion with respect to both the metric and the connections 
the corresponding field equations are obtained as: 



f(R)R^ - \f{R)g^ 



-gf{R) 9 n = 0, 



(16) 



(17) 



in which we have considered the matter action to be in- 
dependent of the connections. From equation (fT7j) we can 
see that the connections are the Christoffel symbols of 
the new metric /i M „ where it is conformally related to the 
original one via the equation 



hftu = f'(R)9iiu- 



(18) 



early universe the equation of state of cosmic fluid is p = 
1/3/9 and the right hand side of equation is zero. The 
Ricci scalar for the FRW metric is as follows: 



on P 2 QR^RR 2 



R 2 



R 2 
n {) 



R 4 



3R 2 Rl + 2i? 4 = 0. (14) 



We divide this equation by Rq and rename y — R/Ro and 
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define a dimensionless variable as r = tR Q . Equation 
(I14[) with the new variables can be written as: 



3y" 



9y' 2 y 



y 



l 



+ y 4 - 3y 2 + 2 = 



(15) 



We expect to have Einstein-Hilbert action for the early 
times of the universe which means y » 1. The nu- 
merical solution of this equation is shown in Fig. [TJ 
Here the Ricci scalar should be larger than Rq but com- 
parable with this term. This assumption implies that 
Ricci scalar should be smaller than H 2 at that time. 
From the FRW metric Ricci scalar is related to the Hub- 
ble parameter through: R = 6H + 12H 2 . Consider- 
ing R(t) as a small and constant perturbation due to 
the modified gravity, the differential equation for the 
scale factor is a + a 2 = ea 2 . By a change of vari- 
able as x = a 2 and considering the initial condition of 
a = at t = 0, the solution of the equation yields: 
x(t) = A sinh( v / e<). Expanding this expression around 
t = results in a(t) = Ct 1 / 2 + at 5 / 2 , where a is a small 
value. 

This deviation of the scale factor from a oc t 1 / 2 depends 
on the Ricci scalar at the radiation epoch. Since R and 
Ro are in the order of Hq 2 , we expect that this model 



Equation (|16p shows that in contrast to the metric vari- 
ation approach (see equation [2]), the field equations are 
second order in this formalism. The trace of field equa- 
tions in the Palatini formalism for f(R) = yj R 2 — Rq 
is: 



R 2 



R 2 -Rl = kT. 



(19) 



The vacuum solution of this equation R = \/2Rq, is the 
same as in the metric formalism. Through the confor- 
mal transformation of the metric g^ v to h^, we get the 
corresponding Ricci tensor: 



Roq = -3- + -/ 
a 2 



a ^ v— 2/q fi\2 ^ , f _i 



(doff - ^VoVo/', (20) 



R tJ = [ah + 2a 2 + f'- 1 {*,} d f + y/'^VoVo/'fe, 

(21) 

where } and V are associated with g M „. Using eqs. 
(fT6|) , ([20j) and (|2T|) we can derive the modified Friedmann 
equation: 

6H 2 + GHf'-'dof + \f- 2 {dQff = <P±M+1 . 

. W 

Now we want to obtain the dynamics of universe in the 
radiation dominant epoch. For this epoch implying the 
equation of state p = 1 /3p, from the equation (fT9|) T = 0, 
we get a constant value of R r = \/2Rq for the Ricci scalar. 
Substituting it in ([22]) results in: 



6H 2 = 2nPr + R ° 



V2 



(23) 
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FIG. 2: Numerical solution of Equations ([25} and (|26|l . Here 
we omit the Ricci scalar between the two equations and ob- 
tain a direct relation between the Hubble parameter and scale 
factor in the logarithmic scale. 



FIG. 3: Numerical solution of scale factor in terms of cosmic 



time, normalized to the R, 



'1/2 



III. EQUIVALENCE OF THE MODIFIED 
GRAVITY WITH SCALAR TENSOR THEORIES 



in which the radiation density changes by the scale factor 
of the universe as p r — poa ~ 4 . The solution of equation 
(1231) results: 



(24) 



This solution is similar to the case of metric formalism, 
showing that for the action of f(R) — \J R 2 — R\ , we 
have almost the same dynamics for the early universe in 
both Palatini and metric formalism. 



For the matter dominant epoch we use the direct de- 
pendence of the Hubble parameter and scale factor on 
f(R) asQ: 



Rf - 3/ 



6/' (-, _ 3 f"(Rf'-2f) 

2 rw-n 



2 i 



(Rf - V) 



(25) 
(26) 



The numerical solution of these equations for f(R) = 
\J R 2 — Rq gravity is shown in Fig. [2l Here we obtain 
the Hubble parameter normalized to its present value in 
terms of the scale factor, in logarithmic scale. By the 
numerical solution of H=H(a) a direct relation between 



the scale factor and cosmic time normalized to the R, 
is shown in Fig. [3J 
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Here in this section we discuss the equivalence of the 
modified gravity with the scalar tensor theories, such as 
scalar-tensor coupled to the curvature like Brans-Dickc 
and non-minimally coupled one which is a result of con- 
formal transformation from the Jordan to the Einstein 
frame [36j. Here in this section we study the correspond- 
ing theories of f(R) = \J ' R 2 — R 2 action. 



A. Brans-Dicke theory 

Theories in which scalar fields are coupled directly to 
the curvature, are termed scalar-tensor gravity [13 ■ Such 
theories can be found in the low-energy effective string 
theory which couples a dilation field to the Ricci curva- 
ture tensor. The simplest and the well known one is the 
Brans-Dicke (BD) theory. The BD theory is a general- 
ization of the general relativity with the action of: 



Sbd 



1 

2k 



uj 



(27) 



where the free parameter loq is often called Brans-Dickc 
parameter and for the case of u>q —> oo GR is recovered. 
The field equation that one derives from the action (|27|) 
by varying with respect to the metric and the scalar field 



2^VVVa</»)(28) 



■5 



V 

24> 



2u>q 



V„V Q c 



R - ^V a 0V a 



<j> - V = Q 



(29) 



where G ^ v = 



T — 

trace of Eq. 



— \Rg^v is the Einstein tensor and 
,ls f is the stress-energy tensor. We take the 
and combine it with Eq. (j2"9")l to omit R, 



which results in the differential equation for the dynamics 
of scalar field as: 



(2w + 3)V a V> = kT + <j)V' - 2V 



(30) 



We can also find an equivalent action for the BD theory 
with introducing a modified gravity as f(R). Let us use 
an auxiliary field \ an d write an equivalent action as: 



S 



1 

2k 



*xy/=g(J(x) + f'(x)(R ~ x)) + S M {9^). 

(31) 

Varying with respect to \ leads to the equation % = R 
if fix) 7^ 0. By redefining the field x> $ = fix) an d 
setting: 



(32) 



= *(*)*- /(*(*)) 
the action takes the form of 

S= hf dix ^~9i®R - Vm + S M (g^,Tp) (33) 



Comparison with the action in ([27)1 reveals that we have 
an equivalent Brans-Dicke theory with loq = 0. In other 
words f(R) theories are fully equivalent to a class of 
Brans-Dicke theories with vanishing kinetic term. For 
the case of our action, f(R) = \J R 2 — Rq we can eas- 
ily find the proper Brans-Dicke potential using equation 



V($) = iW$ 2 - 1 



(34) 



The Ricci scalar also depends on the scalar field by (see 
Fig®: 



R = V' = Rq 



(35) 



For the early universe with large Ricci scalar we have <fi — > 
1 and for the accelerating epoch of the universe, <f> — ► oo 
which results in R — * Rq. Substituting the potential in 
the Brans-Dicke action, the final form of the equivalent 
action with the modified gravity is: 



5' 



BD 



1 

2k 



d 4 x^g-($R - R Vl> 2 - 1) + S M {g^, 

(36) 




FIG. 4: The dependence of the scalar field on the Ricci scalar 
in Brans-Dicke correspondence to the modified gravity. 



Non-minimally coupled scalar-tensor gravity: 
conformal transformation 



Metric theories of gravitation which depend on the 
scalar curvature in a nonlinear way, are usually called 
improperly nonlinear gravity (NLG) models [38| |. A suit- 
able conformal transformation on metric can change the 
NLG lagrangian from the original frame so-called Jordan 
frame into Einstein-Hilbcrt one, minimally coupled to a 
scalar field. It is therefore claimed that any NLG theory 
is mathematically equivalent to General Relativity (with 
the scalar field), but from the physical point of view there 
is a debate in literature if there is an experiment to dis- 
tinguish them [H| . 

To show the equivalence of these two frames one can 
use the auxiliary fields A and B and rewrite the action 
CO) as: 



S 



1 

2i< 



d 4 x^g~(B(R- A) + f(A)) (37) 



Varying the action with respect to B we get R = A which 
results in generalized gravity in Jordan frame. Also vary- 
ing the action with respect to A we get B = f'(A) and 
we can write the action as: 



S=±- 

2k 



d 4 xV=g(f'(A)(R-A) + f(A)) 



(38) 



Now we use the conformal transformation g^ — > 
exp(ip)g tJll , where tp = — In f'(A) and obtain the Einstein 
frame action (scalar-tensor gravity), as follows: 



S = 



K 

A 



d 4 x. 



-g[R 



3 J"(A) 
2 [ f(A) 



m 

f(A)> ^ [f'(A)Y 



(39) 
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FIG. 5: Dependence of potential on the scalar field results 
from the conformal transformation from Jordan to the Ein- 



stein frame in f(R) = \J B? 



R\ gravity. 



We can also write this action as: 

3 
2 



5=4/ <£x^—g[R - \g^d^d vV - V{<p)] (40) 



where the scalar potential is: 

A f(A) 



V{<p) = 



f'{A) f'(Ay 



(41) 



For the case of f(R) = y/R 2 
minimally coupled theory is: 

6 = In ■ 



U 2 



R 



the potential of non- 



VR 2 



R 2 



V(<j)) = Roe-^y/l-e- 2 *. 



(42) 



where for convenience we changed tp to —ip in our calcu- 
lation. For the early universe where R — ► oo, 4> — ► 
and for R — Rq which is the vacuum solution of the 
generalized gravity, <f> = lny/2. Figure (J5|) shows the 
dependence of the potential to the scalar field and the 
evolution of field starts from <f> = to its maximum value 
at <j> = In y/2. The final stage of the scalar field results in 
de Sitter phase expansion of the universe. 



IV. INSTABILITY OF f(R) GRAVITY 

In order to examine the instability of the action we do 
a small perturbation of Ricci scalar around the vacuum 
solution to see if the perturbation grows or decays. A 
small perturbation around the vacuum can be written: 



where R± <C R v and R v is the Ricci scalar of the vac- 
uum. To find the dynamics of R\ , we expand the trace 
of Equation |2]) around R v as: 

f'(R v +R 1 )(R v +R 1 )-2f{R v +R 1 )+3WV IM f'{R v +Ri) = 0, 



where the first order perturbation gives: 
Rvf"{Rv) — f (R 



Ri = 0. 



For the case of f(R) = JR 2 



R 2 , we have: 



V^V„i?i + bR x = 0, 



(44) 



(45) 



(46) 



where b — ^j2i? . For FRW metric we can write this 
equation as follows: 



Ri 



\t){dt 



8%)Ri + bRx = 0. (47) 



Since the vacuum solution of the background results in a 
de Sitter solution and exponential expansion of the uni- 
verse, we use a(t) = exp2 a (' -t °) for the dynamics of the 
background, in which a = ( ^ 2 ^-° ) \ . In the Fourier space 
(1471) can be written as: 



Ri(k) + tfe-^-^Rtfr) - bRi(k) = 0. 



(48) 



It is clear that in the limit t — > oo we can ignore the 
second term and the solution is exponentially growing. 
We can also do semi-analytical calculation by using a ~ 
t^ 1 and b ~ t^ 2 and as a result, the differential equation 
is: 



R"i(k) + {e x ~ T k 2 - l)i?i(fc) = 0, 



(49) 



where prime represents derivation with respect to the 
dimensionless variable r = ^ . For perturbations smaller 
than horizon, (i.e. k < 1) we will have exponentially 
growing modes, while wavelengths larger than horizon 
will oscillate. Figure © shows the dynamics of R\ in 
terms of the normalized time r for the modes smaller 
than horizon. 

In order to stabilize this action we can add a quadratic 
term to the the action as aR 2 and examine the stabil- 
ity. This stability condition can be also obtained if one 
directly solves the perturbed generalized Einstein equa- 
tion: 

- R'l + e- (T - 1} (^ + d 2 y + d 2 z )R l + bRi = 0, (50) 

where b = ^{Rq 1 ~ 2a) _1 . For r < 1 and in the 
Fourier space, one gets 

Ri(k) + (ek 2 - b)Ri(k) = 0. 



(51) 



So Ri(k) would oscillate if ek 2 > b, or equivalently, 
the model would be stable for perturbations of all wave- 
lengths if 



R = Ry + i?l , 



(43) 



2R a > 1 



(52) 
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FIG. 6: Numerical solution of Equation (|49[l for the pertur- 
bation of Ricci scalar around vacuum solution. 



Adding aR 2 term not only makes the action stable but 
it will be also important for the early time inflation. To 
show this we can ignore the first term of action at the 
early universe while R » Rq and consider the action 
as: 



f(R) ~ aR 2 



(53) 



For the radiation dominant epoch where p = 1/3/9, the 
trace of energy momentum tensor is zero, = 0. So 
the generalized Einstein equation reduces to 



VVJ'(fl) = 0. 



(54) 



For homogenous universe we will have only temporal 
derivative of f'(R) — 2aR, which results in: 



V°V i? = 0, 



(55) 



where for FRW metric Tq A = and we can replace the 
covariant derivatives with partial ones which results in: 

R = Ct + D, (56) 

C and D can be chosen so that when t ~ 0, R ~ D, 
So the scale factor would behave exponentially (a(t) oc 

e"^*) in the very early epoch. 

The same result can be achieved by looking to the in- 
stability of scalar potential of V{4>) in conformal trans- 
formation of metric to Einstein frame. For the case of 
f(R) = \/ R 2 — Rq the corresponding scalar potential as 
shown in Figure ([5]) is instable at the maximum of the 
potential since V < 0. Adding the quadratic term gives 
the scalar potential as: 



V(4>(R)) = 



R 



y/R 2 



R 2 



aR 2 



2aR ( 



2aRf 



(57) 



Using the trace of generalized Einstein gravity, vacuum 
happens again at R = V2Rq. The second derivation of 
potential at R — \/2Rq yields: 



V"(R= V2R ) = 



-1 



AR 2 a 2 



R {l + 2R a) 4 



(58) 



The stability condition V" > 0, implies 2Roa > 1, which 
is in agreement with (|52[) , the approach of perturbation 
of Ricci scalar. 



V. CONCLUSION 

Here in this work w e introdu ced a new action for the 
gravity as f(R) = \J R 2 — R\ which can imply a late 
time acceleration for the universe. Expanding this action 
results in F(R) = R — 1/2Rq/R + ... which is similar to 
that of 1 1 R gravity models. In the strong gravity regime, 
action can be written as f(R) = R and the Einstein- 
Hilbert action is recovered. By the two different metric 
and Palatini approaches we obtained the field equations 
of gravity and showed that by tuning the parameter of 
model Rq , universe can enter an acceleration phase in 
any desired redshift. 

We also obtained the equivalent theories as the Brans- 
Dicke and non-minimally coupled scalar-tensor gravity. 
Finally the instability of this action was examined by 
perturbing the Ricci scalar around the vacuum solution 
which showed that we have instability for the perturba- 
tions smaller than the horizon wavelength. Adding the 
quadratic term of aR 2 to the action stabilize it for this 
type of perturbations for 2aRo > 1. We also tested the 
stability condition for the scalar field in the equivalent 
scalar-tensor gravity which resulted in 2aRo > 1 . 



Appendix 

In Section [TT| the solution of the generalized Einstein 
equation in the spherically symmetric space for the spe- 
cial case of vacuum solution is obtained. We try to find 
another solution perturbing Ricci scalar around the vac- 
uum solution as: 



R = V2R a + f(r)e, 



(59) 



where e is a small dimensionless parameter. Substituting 
this term in the vacuum field equation of ((4} and ignoring 
the terms of second order in £ gives: 



f(r)" - ^Rof(r) = 0, 



(60) 



where the solution is as: 



f(r) = Asin(\ 



'2V2Ro 



Bcos (\ 



2V2Ro 



r). (61) 
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We can ignore the "sin" term by the argument that 
the maximum of Ricci scalar is around the gravitational 
source which should be vanished at infinity. So the over- 
all Ricci scalar can be written as: 



(62) 



Substituting the Ricci scalar in ([5]), we obtain the metric 
as: 



A(r) = - ln[l + 



Cl c 2 \/2 2 

2 ~ To" ° r 

r r z 12 



(63) 



If 1 3 

-Roei^ / r 3 cos (-2* s/m r)dr 



It is obvious that if we set e = (ignore the perturbation 
term) we will achieve the solution of constant Ricci scalar 
and we can determine c\ = —2GM/c 2 and c-i = . 
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